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Production rates of dileptons and photons from the quark-glnon (QGP) phase are calculated 
taking into account effects of confinement and spontaneous chiral symmetry breaking (xSB) not 
far from the transition temperature Tc- We find that the production rates of dileptons with large 
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temperature range and the steepness of the crossover, are reflected in these quantities. 
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I. INTRODUCTION 

In the limit of massless quarks with Nf flavours, 
quantum chromodynamics (QCD) has an exact chiral 
SU{Nf)p X SU{Nf)ii symmetry. Dynamical quark mass 
generation implies that this symmetry is spontaneously 
broken. In addition, non-zero current quark masses in 
the QCD Lagrangian break chiral symmetry explicitly. 
Spontaneous chiral symmetry breaking (ySB) is manifest 
in the temperature dependence of the chiral (quark) con¬ 
densate, {qq). Recent determinations of {qq)T from lat¬ 
tice QCD with Nf = 2-1-1 flavours and physical current 
quark masses display a chiral crossover at a transition 
temperature around ~ 155 MeV^ [1, 2]. Schematic 
models of the Nambu & Jona-Lasino (NJL) type [4] link 
the dynamically generated quark masses, rn{T), to the 
chiral condensate: 

MT) {qq)T ,, .. 

m{0) ^ {qq)o • ^ ^ 

This concept is adopted in the present work, as will be 
explained in the next section. 

The other important aspect of QCD at temperatures 
around and above the transition temperature (Tc) is the 
deconfinement transition. At zero temperature quarks 
and gluons are confined inside hadrons. Above the de¬ 
confinement temperature the quarks and gluons are re¬ 
leased and become active degrees of freedom. An order 
parameter for this transition to the quark-gluon plasma 
(QGP) is the Polyakov loop. In the pure gauge case 
(Nf = 0) deconfinement emerges as a first-order phase 
transition, whereas it becomes a continuous crossover in 
the presence of dynamical quarks. Thus, in full QCD, 
the Polyakov loop is not an order parameter in the strict 
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^ This updated value of Tc, slightly lowered from a previous de¬ 
termination [3], resulted from improvements in controlling the 
continuum limit. 


sense but it nonetheless serves as a monitor for decon¬ 
finement as a rapid crossover transition. Motivated by 
this scenario, models have been designed (referred to as 
PNJL models) [5] that unify the Polyakov-loop-induced 
suppression of color non-singlet degrees of freedom be¬ 
low Tc with the NJL mechanism of spontaneous chiral 
symmetry breaking. 

High-energy heavy-ion collisions (HIC) are the experi¬ 
mental tool for investigating properties of the QGP. Elec¬ 
tromagnetic probes such as dileptons and photons [6-13] 
produced in these collisions are particularly important 
messengers for the physics of the hot and dense QCD 
matter since they are barely affected by successive inter¬ 
actions on their way through the surrounding expanding 
medium. At the Relativistic Heavy Ion Collider (RHIC), 
the PHENIX experiment has measured [10] a large ellip¬ 
tic anisotropy of the produced photons characterized by 
the quantity V 2 - Since the anisotropy is considered to be 
produced by the collective flow at later stages of the HIC, 
it is natural that the produced hadrons, mainly generated 
at the later stages, have large V 2 - By contrast, photons 
are generated also at early stages of the HIC (in the ini¬ 
tial state [14], or in the thermalization process [15]). It is 
therefore not easy to explain [16] the large photon V 2 , and 
currently no theory has completely succeeded in doing it. 
Recently, a new mechanism to generate a large photon V 2 
was suggested [9]. It is based on the suppression of the 
photon production in the QGP phase due to the confine¬ 
ment effect. Nevertheless, the ySB effect was not taken 
into account in this analysis, so if this effect further sup¬ 
presses the photon production from the QGP phase, one 
can expect that the photon V 2 is more enhanced. 

Motivated by these issues, we calculate the production 
rates of dileptons and photons in the QGP phase, using 
a model in which the confinement and the ySB effects 
are both taken into account. The present work focuses 
in this context primarily on the role of the latter effect, in 
terms of dynamically generated constituent quark masses 
around Tc. It is an extension of the previous Refs. [9, 17] 
in which corresponding studies have been performed with 
quark masses set equal to zero. 
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This paper is organized as follows: The next section 
focuses on the implementation of the effects of confine¬ 
ment and in our model. The dilepton production 
rate is calculated and the xSB effect on this rate is dis¬ 
cussed in Sec. III. Section IV is devoted to the calcula¬ 
tion of the photon production rate, and to the discussion 
on how the xSB affects this rate. For this purpose, the 
photon production rate in the case of massless quark is 
also calculated beyond leading-log order^. We summarize 
and present concluding remarks in Sec. V. The detailed 
derivation of the photon production rate for vanishing 
quark masses, m = 0, is given in the Appendix. 


II. IMPLEMENTING CONFINEMENT AND 
CHIRAL SYMMETRY BREAKING EFFECTS 

This section describes the treatment of confinement 
and spontaneous chiral symmetry breaking within the 
framework of the present model. 


A. Polyakov loop 


The confinement effect is taken into account as a modi¬ 
fication of the quark and gluon thermal distribution func¬ 
tions, in a so-called semi-QGP model [9, 17-20]. The 
starting point is QCD with Nc colors and with finite aver¬ 
ages of the temporal component of the background gluon 
field (A°), given by: 


A° = -diagiQ\Q^...,Q^^), ( 2 . 1 ) 

9 


where g is the coupling. A° is assumed to be diagonal in 
color space and constant in space and time. The spatial 
components of the background gluon field. A®, are set 
equal to zero. The tracelessness of the gauge held implies 
(5“ = 0. Neglecting fluctuations of Aq, this quantity 
is related to the expectation value of the Polyakov loop 
(Oby 


Z=—TrL=—y, 


,i/3Q“ 


( 2 . 2 ) 


where 


L(x) = ^Pexp dTgA°(T,x)^^ (2.3) 


is the Wilson line in the temporal direction in imaginary 
time, with /3 = l/T, and V denotes path ordering. 

When performing calculations involving quarks and 
gluons, the background gluon held acts as an imaginary 
chemical potential coupled to color charges. The quark 
and gluon distribution functions, np{k^) = (e^^ -I- 1)“^ 
and nB{k°) = — 1)“^, are modihed as [19] 

npaik ) = ^) (2-4) 

nBab{k ) = (2-5) 

where a and b are color indices in the double-line nota¬ 
tion [19, 21], running from 1 to and <5“^ = Q°' — Q^- It 
might appear that the distribution functions have imagi¬ 
nary parts, but in fact, after summing over color indices, 
these distributions turn out to be real as they should. To 
see this, consider the color-averaged quark distribution 
function: 


-| OO 

(npaik^)) ^—J2nMk°) = , 

^ a n—1 

( 2 . 6 ) 


where we have introduced In = TtL'^/Nc- Since In is real, 
it follows that {npaik'^)) is real. 

Next we demonstrate the suppression of the distribu¬ 
tion functions for color-non-singlet degrees of freedom in 
the conhned phase. From Eq. (2.6), the averaged quark 
distribution function becomes 

{nFa{k°)) = _|_ I ’ 

where we have used the following expression, valid in the 
confined phase for odd^ N^'- 


I (_i)fc(Afc-ei) (^n = kNe) 

\ 0 (otherwise) . 


( 2 . 8 ) 


Non-trivial contributions to in the expansion (2.6) 
arise whenever n is an integer multiple of (odd) Nc- We 
see that (n_Fa(fc°)) is suppressed compared to the distri¬ 
bution in the deconfined limit, = 0, and it vanishes 
at large Nc- For Ac = 3 one notes that {npaik'^)) repre¬ 
sents color-singlet 3-quark objects (with energy E = 3k^ 
and baryon number i? = 1), although these structures 
are not localised or clustered in space. Sometimes this 
property is called “statistical confinement”. The gluon 
distribution function is likewise suppressed after taking 
the color sum. 


^ Leading-log approximation means: regarding a quantity of order 
one as subleading compared with a quantity of order ln(l/e) if 
e is very small. In the case of photon production, e = g'^TjE, 
where E is the photon energy, g the coupling constant, and T the 
temperature, respectively. The production rate in the massless 
case at leading-log order was already calculated in Refs. [9, 17]. 


® The expression remains valid also in the case of even Nc, but with 
an unphysical result. For example, when Nc = 2, the baryonic 
excitations appearing in the confined phase are diquarks, i.e. 
bosons, whereas the distribution function Eq. (2.7) is fermionic. 
For this reason, we restrict ourselves to the odd-Ac case. 
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FIG. 1: Left panel: The Polyakov loop (1) deduced from lattice QCD computations [3], and the reduced Polyakov loop in 
which the perturbative correction is removed as a function of T. Right panel: temperature dependence of Q. Here 

^MS' ~ Tc/1.35, where Tc = 155 MeV [1, 2]. The curves In both panels connecting the points are drawn to guide the eye. 


For the Nc = 3 case can be written in the repre¬ 
sentation 

(Q“) = (-g,0,Q), (2.9) 

where Q with 0 < Q < 2ttT is the only independent 
quantity. The Polyakov loop reads: 

i = ‘ + (2.10) 

In the asymptotic deconhned phase {I = 1), we have Q = 
0, while in the confined phase, Q = 27rT/3 for the pure 
glue case (1 = 0). A correspondingly similar expression 
holds for Z„: 


1-1-2 cos{nPQ) 


( 2 . 11 ) 


In the present work, and Q are determined from lat¬ 
tice QCD data in the following way [20]: First, we remove 
the perturbative correction [22, 23] from the Polyakov 
loop, focusing on nonperturbative effects: 


apparent convergence” criteria [23, 24]: 


= 247r^ 




+ Nf (A\n2- 1) - 


llNr 


1 -1 


(2.14) 


mjj = {2Nc + Nf)ATr^T^ (llNc - 2Nf) 

.O- 


lE 
2 I At'2 


-h4Al/ln2- 


bNl +N]+ 0Nf/(2N,) 


2N, + Nf 


-1 


(2.15) 


where is the renormalization scale in the modified 
minimal subtraction scheme and — 0.577 is Euler’s 
constant. From Eq. (2.12), we see that 1{Q = 0) exceeds 
unity due to the renormalization effect. Now we assume 
the following relation which reduces to Eq. (2.12) at Q = 
0 and 51 1: 


l(Q) = e^^^^=°U^°\Q). (2.16) 


1{Q = 0) 
6l(Q = 0) 


1 -l- 6l(Q — 0), 


8ttT 


9*Cf 

(47r)^ 



+ N, 



( 2 . 12 ) 


(2.13) 


where Cf = (N^ — l)/(2Ni.) and mp is the Debye mass 
of the gluon. We use the running coupling constant and 
the expression of the Debye mass estimated by “fastest 


With I deduced from lattice QCD [3], we can determine 
from Eq. (2.16), and then obtain Q from Eq. (2.10) 
using . These quantities are plotted in Fig. 1, set¬ 
ting = Tc/1.35. Note that is still different from 
unity even around ^ 3Tc, where Tc is the pseudo-critical 
temperature, approximately 155 MeV in recent lattice 
computations [1, 2]. 

We add the following remark concerning the interpo¬ 
lating ansatz, Eq. (2.16): If the perturbative correction 
is not considered (i.e. SI = 0), then 1{Q) gives the up¬ 
per bound for the Polyakov loop. On the other hand, 


















4 


using Eq. (2.16) for not too small 61 presumably overesti¬ 
mates the perturbative correction and is expected to give 
a lower bound for . 


B. Dynamical (constituent) quark mass 


TABLE I: Parameters (in MeV) used in the expression (2.17) 
for the dynamical quark mass m(T). 


Set 

A 

B 

C 

D 

S 

22 

22 

13 

33 

To 

180 

150 

180 

180 


Quark mass generation in QCD implies spontaneous 
ySB. The dynamical (constituent) quark mass m is not 
directly accessible in lattice QCD computations, but it 
figures as a well-defined quantity in studies of the quark 
propagator using Dyson-Schwinger equations [25] and 
schematic models of NJL type [4]. Starting from massless 
current quarks, such approaches feature, in the form of 
a characteristic gap equation, a proportionality between 
the dynamical quark mass m and the chiral condensate 
{qq), realised also at finite temperatures T as indicated 
in Eq. (1.1). The transition from the ySB (Nambu- 
Goldstone) phase to chiral symmetry restoration in the 
Wigner-Weyl phase with (qq) = 0 is second order in the 
limit of vanishing current quark masses. This transition 
becomes a continuous but rapid crossover when the ex¬ 
plicit chiral symmetry breaking by small non-zero current 
quark masses is included. 

In the present work we make use of the typical be¬ 
haviour of m(T) derived from NJL and PNJL models. 
The input in such models is designed to reproduce vac¬ 
uum properties of the pion (its decay constant and its 
physical mass starting from non-zero current quark 
masses). In the two-flavour case the emerging dynami¬ 
cal u and d quark masses can be well represented by the 
following expression: 


m(T) 


mp 

2 


1 -arctan 

TT 


(T-To 

\ s 


(2.17) 


A standard fit to two-flavour NJL results gives mp ~ 350 
MeV, J ~ 20 MeV and a transition temperature Tp ~ 180 
MeV. The PNJL model features a characteristic entangle¬ 
ment of m(T) reflecting spontaneous (dynamical) ySB, 
with the Polyakov loop representing confinement. While 
values of Tp found in PNJL calculations are typically 
larger than those found in NJL models, the correspond¬ 
ing profiles of m(T) are qualitatively similar in all such 
approaches, apart from shifting the T-scale [5, 33]. 

The Nf = 2 NJL and PNJL model calculations yield 
values of Tp that are usually higher than the transition 
temperature determined in recent Nf = 2+1 lattice QCD 
computations [ 1 , 2]. Since one of the aims of the present 
paper is to explore the systematics of quark mass effects 
on dilepton and photon production rates, we allow varia¬ 
tions of Tp between 150 and 180 MeV, and we also exam¬ 
ine the influence of different slopes of the chiral crossover 
controlled by the parameter 6. Pour parameter sets will 
be used, as given in Table 1. The constituent quark mass 
Too at T = 0 will be kept fixed (at mp = 347 MeV) for 
all sets A - D. Set B with Tp = 150 MeV is motivated by 
the chiral transition temperature ~ 155 MeV found 
in the lattice QCD computations of Refs. [1, 2]. Sets 



T [MeV] 


FIG. 2: The dynamical (constituent) quark mass, Eq. (2.17) 
as function of T, evaluated using parameter sets A, B, C, and 
D. Also plotted for comparison is the thermal quark mass, 
Eq. (2.18). 


C and D with different values for 6 are used in order to 
analyze the sensitivity to the steepness of the crossover. 
With these parameter options the behavior of rn(T) in 
both NJL and PNJL models can be covered, their main 
difference being a shift in Tp. 

The temperature dependent quark masses m{T) with 
these parameter sets are shown in Fig. 2. For later con¬ 
venience, we also plot the thermal quark mass in the 
deconfined phase (with W = 3): 


mth{T) 


iL 

ye ’ 


(2.18) 


with g determined by Eq. (2.14). The dynamically gen¬ 
erated quark mass tends to exceed the thermal mass at 
temperature T < 200MeV. 


III. DILEPTON PRODUCTION 

At leading order of the strong and electromagnetic cou¬ 
pling, the dominant process contributing to the dilepton 
production rate is the quark-antiquark pair annihilation 
process followed by pair production of the leptons via the 
photon. The amplitude of this process factorizes into a 
part describing the gg —> 7 process (shown in Fig. 3), 
and a part describing the conversion of the photon into 












(3.4) 



FIG. 3: Quark-antiquark pair annihilation process contribut¬ 
ing to dilepton production. The solid lines represent the quark 
or antiquark while the wavy line represents the (timelike) pho¬ 
ton that converts into the lepton pair. 


the dilepton [26]. The differential production rate is'*’ 

S (3-1) 

where pi and p2 are the four-momenta of the leptons, p = 
P1+P2, and a = e^/(47r). Here n‘'f^^{p) is the Wightman 
photon self-energy, which includes the information about 
the dynamics of the gg —1 7 process. The leptons are on- 
shell {pI = P 2 = 0 neglecting the small lepton masses), 
hence p is time-like {p^ >0). At leading order we have 
the following expression for 7T<^^(p): 




- s / 

a,/,spin 


2Ei{2nf 


2E2{2TTf 




X {2'k)‘^5^‘^\p - ki- k2) na{Ei)na{E2), 

(3.2) 


where / is a flavour index running from 1 to and 
k^ = Ei = i/k^ -I- m?{T) are the quark and antiquark 
energies, m(T) is the temperature-dependent constituent 
quark mass and na{E) = [ed(^+®Q“) -p l]“i is the distri¬ 
bution function for anti-quarks [19] with a sign change of 
the imaginary chemical potential as compared to 
na{E) = [ed(^“®Q ) -p l]“i for quarks. The squared ma¬ 
trix element is 

^ = 4e2q'^(p^-b 2m^), ^3 3^ 

spin 


where we have used fcj = fc| = m^, and qf are the electric 
quark charges in units of e for each flavour /. 


A. The jpj = 0 case 

Consider first the case with zero total three-momentum 
of the lepton pair. When jpj = 0 the integration in 
Eq. (3.2) can easily be performed. The leptons are pro¬ 
duced back-to-back. The energies of the quark and anti¬ 
quark satisfy Ei = E 2 = E = p°/2. As a result Eqs. (3.1) 


Note that the differential rate dF is defined so that it is propor¬ 
tional to the phase space element (d^pi/2Ei) (d^p2l2E2) of the 
lepton pair, and therefore carries mass dimension four. 


dr 

d'^p 


127r4A2 




y/E'^ — w? 
E 



X ^ na{E) na{E) 0{E - m). 

a 


In the case E < m the dilepton production rate vanishes 
for obvious kinematic reasons. The physical interpreta¬ 
tion of this expression is evident. It can conveniently be 
factorized into three parts, coming from the phase space 
integral, the matrix element squared, and the distribu¬ 
tion functions. For Nc = 3 Eq. (3.4) becomes 


X 9{E — m) , 

where 


A{m) = 
B{m) = 


C{Q) = 


^E^-m?{T) 

E 

mP'lT) 

1 -I- —- 

2E^ ’ 

e-2/3£ + [;(g) 


1 ]' 


e-2/3£ + [3^(Q) _ 1] e 


■/3fS_gl 
-PE + l 


come from these three parts, respectively. Note that A 
and B become unity when chiral symmetry is completely 
restored (m = 0), while (7=1 when there is no confine¬ 
ment effect {I = 1). For m = 0 our expression reduces to 
that in one of the author’s previous works [9, 17]. 

Only the parts A and B depend explicitly on the quark 
mass m[T). The term C depends on the quark and an- 
tiquark energy E, but this is in turn fixed by the lepton 
energy, i.e. E = p°/2. We also note that the product 
AB is always less than one as long as > 2m. This im¬ 
plies that finite m always reduces the dilepton production 
rate. By contrast, the background gluon field Qa enters 
directly in part C through the sum over products of quark 
and antiquark distributions, na{E)na{E). This separa¬ 
tion of terms in Eq. (3.5) helps making the analysis of 
ySB and confinement effects transparent®, especially in 
comparison with previous results of a semi-QGP model 


® Here the dependence on the number of flavour enters only 
through the prefactor 57 / 9/1 because we have assumed a com¬ 
mon T-dependent dynamical quark mass m for all (light) flavours 
throughout the present work. For Nj = 2 we have 57/9/ = 3/9. 

® Given the entanglement of confinement and xSB as it is realized 
in the PNJL model, m and Q are not independent quantities. It 
is nonetheless instructive to analyze the effects of confinement 
and xSB on the dilepton production separately, and so it is con¬ 
venient in the present context to treat m and Q as independent 
quantities. 
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FIG. 4: The quark mass dependence (displayed by the quantity AB) and the Polyakov loop dependence (displayed by C) of 
the dilepton production rate at |p| =0, as functions of temperature T. Also shown is the T dependence of the product ABC. 
Parameter set A is used to obtain the input mass m(T). In the left panel, is set to 700 MeV while in the right panel, = 1 
GeV. The curves in both panels connecting the computed points are shown to guide the eye. 



FIG. 5: Results for the quantity AB using the four parameter 
sets (A, B, G, and D) of Table I. The energy is chosen to 
be 700 MeV. The curves connecting the computed points are 
drawn to guide the eye. 


calculation [9, 17] in which the xSB effect was not con¬ 
sidered. The difference is then primarily in the function 
AB. 

The quantities AB and C are shown as functions of T 
in the left panel of Fig. 4. In the given example is fixed 
at 700 MeV, where we have used parameter set A for the 
evaluation of the quark mass m. The AB plot displays 
the m dependence while the latter plot contains informa¬ 
tion about the dependence on 1. From the AB plot one 
observes that this quantity is significantly smaller than 
unity in the region T < 180 MeV: the constituent quark 
mass suppresses the dilepton production rate by modi¬ 


fying the available phase space and the matrix element. 
When T > 180 MeV AB quickly approaches unity be¬ 
cause around this temperature, m decreases rapidly as 
shown in Fig. 2 and becomes negligible compared to the 
energy E = p^/2 = 350 MeV appearing in AB. From 
the plot of C one observes that the effect of the Polyakov 
loop does not suppress the dilepton production, but even 
slightly enhances it. This unexpected behaviour has al¬ 
ready been discussed previously in Refs. [9, 17], referring 
to the fact that the initial state (the quark-antiquark 
pair) is color-singlet. 

The product ABC indicates the modification of the 
dilepton production rate compared to the case with I = 1 
and m = 0. It shows the combined effect of m and the 
Polyakov loop. Since C > 1, the suppression coming 
from finite m in the region T < 180 MeV is balanced 
by the enhancement due to the Polyakov loop. For T > 
180 MeV the factor AB approaches unity and ABC is 
nearly equal to C. With increasing delepton energy 
the influence of m becomes less significant. This is illus¬ 
trated by comparing the left and right panels of Fig. 4. 
At = 1 GeV, AB differs only slightly from unity even 
when T < 180 MeV. 

We also examine the dependence of the dilepton pro¬ 
duction rate on the four different parameter sets summa¬ 
rized in Table I. The result for AB is shown in Fig. 5 
using = 700 MeV. Lowering Tc increases the dilep¬ 
ton rate as the comparison between “set B” and “set A” 
results shows. This trend is expected since, as Tc is low¬ 
ered, m becomes smaller in the relevant range and so the 
suppression effect on dilepton production is reduced. On 
the other hand, if the chiral crossover transition proceeds 
with a steeper slope around Tc^ the quark mass m below 
the transition temperature tends to be larger and the 
suppression effect induced by m is enhanced, as can be 
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seen from the comparison of the plots for “set A”, “set 
C”, and “set D”. Above the transition temperature the 
quark mass quickly decreases and the effect is not visible 
in the figure. 


B. The IpI A 0 case 


For non-vanishing three-momentum of the lepton pair, 
IpI 7 ^ 0, using 

i5(p° — El — E 2 ) = ■^^-q-fT-^(5(cos0 — cos^o), (3.6) 

|Pl|ki| 


Kinematics implies that the dilepton production rate 
vanishes for < 4m^. Performing the integration one 
arrives at: 


dr 

d'^p 


12^ E e{p^ - 

1 -I- g-PiE+-iQa,) \ 

1 -I- Q-P(E--iQa) J 
g-/3p“ _|_ ^-P{E+-iQo.) 

g-/3p0 _|_ g-p(E_-iQa) 

(3.8) 





where cos^o = —(p^ — 2p°i?i)/(2|p||ki|) and 0 is the 
angle between p and ki, Eqs. (3.1) and (3.2) give: 


dr 

d'^p 


2 

12'k‘^p'^ IpI 


0{p^ 


L 


iS+-p “/2 / 0 

dEUa 

E--p°/2 


■E 


4m^) 




where 


E± = 



In contrast to the |p| = 0 case, the effects of the Polyakov 
loop and of the quark mass m are not well separated any 
more. In the limit m = 0 we reproduce the corresponding 
expression in Ref. [9, 17]. 

In the Nc = 3 case, we have 


dr 

d'^p 


„2 


127r4 2 :^ 
/ 


qj nsip )0ip - 


X 3 


1 _(_ g-/3(p -l-|p|)/2 


T 

^^^M^'"l + e-/3(p“-|pl)/2 


B(rn)ffi{Q,m), 

(3.9) 



^IpI ±p° 




where B{m) = {p^ + 2m'^)/p^ and fij{Q,m) = 

fij{Q,m)/fij(p,0), with 


/„(e,ra) = inlYT3iFwv 


- .“t/ P-2/3S+ 


Note that B = 1 when m = 0, and fij{Q,m) becomes 
unity when I = 1 and m = 0. 

The effect of the Polyakov loop I is analyzed in Refs. [9, 
17]. Here we focus on the role of the dynamical quark 
mass, m(T). The quantity B{m)fij{Q,m)/fij{Q,0), 
which is the ratio of the dilepton production rate at finite 
m to that at m = 0 {dBjd'^pj{dr /d'*p|m=o))i is shown in 
Fig. 6 as a function of |p| at three temperatures around Tq 
(140, 180, 190 MeV), and at one temperature (T = 540 
MeV) that is significantly larger than Tq. In the left 
panel, p° is fixed at IGeV and the parameter set A is 
used. Suppression effects on dilepton production due to 
non-zero quark masses are marginal at low momentum 
IpI but become significant at higher momenta. Complete 
suppression occurs at the limit, = 4m^, due to the 
kinematic constraint. At lower temperature the quark 
mass m{T) is larger than at higher temperature, hence 


-3;e-2/5p“g-/3£^+ 


(3.10) 


the suppression is more significant. 


The parameter dependence of the dilepton production 
rate encoded in B{m)fij{Q, m)//^{Q, 0) is plotted in the 
right panel of Fig. 6 as a function of |p|, at T = 170 MeV 
and p° =1 GeV, using the parameter sets A, B, C, and 
D. The differences in m{T) translate directly to the mo¬ 
mentum (pfc) at which B{m)fij{Q,m)/fij{Q,0) vanishes 
because of the kinematic constraint: Pk\c < Pk\ A < Pkloi 
and Pkls is larger than pk\A- Since pk = \J (p°)^ — 4m^ 
is determined entirely by m{T), this tendency is under¬ 
stood in the same way as explained in the discussion of 
the quantity AB in the previous subsection. 
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p [MeV] 


FIG. 6: Ratio of dilepton production rate, dF/d'^pl{dr/d'^p\jn=o), with and without inclusion of dynamical quark mass m{T), 
as a function of |p| with =1 GeV. Left panel: using parameter set A. Results are plotted at different temperatures. T = 140 
MeV: solid (red) curve; T = 180 MeV: dashed (black) curve; T = 190 MeV: thick solid (blue) curve; and T = 540 MeV: dotted 
(magenta) curve. Right panel: same quantity at fixed temperature T = 170 MeV for different parameter sets. Parameter set 
A: solid (red) curve; set B: dashed (green) curve; set G: thick solid (blue) curve; set D: dotted (magenta) curve. 


IV. PHOTON PRODUCTION 

At leading order the dominant contributions to the 
photon production rate come from Compton scattering 
{qg —>■ q'y) and pair annihilation {qq —^ 57) [9, 17, 27-29]. 
The Landau-Pomeranchuk-Migdal (LPM) effect [30-32] 
also contributes at the same order when m = 0. In this 
paper we neglect the contribution from the LPM mecha¬ 
nism. This is justified when m is large compared to the 
thermal quark mass (m ^ mth ~ gT). It is also justi¬ 
fied, even when m = 0, in the large Nc limit [17]. The 
two basic processes, quark-gluon Compton scattering and 
pair annihilation of quarks, are illustrated in Fig. 7. The 
amplitudes of these processes are infrared-singular in the 
limit of small energy and momentum of the exchanged 
quark. These singularities are naturally regularized when 
the constituent quark mass is finite. The resulting photon 
production rate has a characteristic logarithmic depen¬ 
dence on m as will be shown later. Even when m = 0 
the thermal quark mass acts to regularize the amplitudes 
such that they depend logarithmically on mth(T). Which 
mass acts prominently as a regulator depends on which 
one is larger at a given temperature. With the purpose 
of this paper to investigate finite-mass effects systemati¬ 
cally, we calculate the photon production rates with finite 
dynamical quark masses and compare versus results with 
m = 0 but using a resummed quark propagator, which 
contains the information on the thermal quark mass. 


A. Compton scattering 

Consider first the contribution from the Compton scat¬ 
tering process. Within the kinetic theory, the photon 
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FIG. 7: Basic photon production processes: Compton scat¬ 
tering and pair annihilation. The curly line represents the 
gluon. 


production rate from this mechanism is expressed as 



X 

X 


\M\^ 


nFa{E\^) nBcd{\^'\) [1 - 


k' - p - p') 
(4.1) 


where the sum is over quark flavours /, the spins of all the 
particles, and color indices {a,b,c,d). M is the matrix 
element of the Compton scattering process, the explicit 
expression of which will be given later. We consider a 
real photon with p'^ = \p\ = E and focus on the E ^ T 
case relevant for high-energy heavy-ion collisions. The 
Boltzmann approximation will be used in the initial state: 
in Eq. (4.1), we replace nF{E\^-iQa) nB{\^'\-iQc+iQd) 
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by exp[-/3(£'k + |k'|-iQa-*Qc + *Qd)]- This is justified 
when Ei^, |k'| ^ T. It is known that this treatment gives 
the correct result at leading-log order [28, 29]. With this 
approximation Eq. (4.1) leads to the following expression: 


E 


dr 

d^p 


Comp 


1 T 
32{2TrfE 


e-PE 




^ g^/3(Qa—Qb+Qc—Qd) 

X In -I- g-/5(™^.E/s+s/(4£;)-iQfc)^ 

(4.2) 


Here we have introduced the Mandelstam variables: s = 
{k + k'Y and t = {p — k)^. The variables with tilde are 
defined as s = s — and t = t — m ?, and we have intro¬ 
duced t+ = —m^sls and t- = —s. For details concerning 
this change of integration variables, see Refs. [17, 27]. 

The squared matrix element after spin summations is 
given by: 


J2\Mf = 32e^g\}[{E%a? 


spin 


2 ~ r\ ^ 

- \- m s + 2m 


+4^ f “ w+ 

V 2 


st 


[s + i + 4m^) 

(4.3)" 


where {r<^)ba = with 

encodes the color structure of the quark-gluon vertex in 
double-line notation. The first term is the contribution 
from the s-channel process, the second from t-channel 
exchange, and the third term comes from the interference 
of these two amplitudes, respectively. One finds: 


E 


dr 


Comp ^ f 

nOO 

X / dsG{s,m^) 

Jo 

1 p-nP{ndE/s+s/{iE)) 

*'71. C' 


E 


(_i)"+i 


where Us = J (47r) and 


G(s, rr?) = 


s r 2m} _ vri} 

4(s -I- m^Y' a 

, 1 rr? Am}' 

- -- + 2^ + — 


)ln 

s 

1 H- 2 

/ 



(4.4) 


(4.5) 


In the deconfined phase with /„ = 1, this expression be¬ 
comes: 


E 


dr 


X G(s, m^) In (l -h , 


ds 


(4.6) 


In the limit m} —> 0 the logarithmic term in G(s,w}) 
dominates over the other pieces: 


^ djp 


Comp 


IGtt^ 




T 


-PE 





ds 


xE 


(zlZEg-n/3s7(4B) 


1 


S 

m? 




f 


^ (-1)^+^ , , FT 

n—l 


(4.7) 


This expression agrees with the one in leading-log or¬ 
der when regarding as the infrared cutoff /i^ used in 
Ref. [17]. It also suggests that m works as an infrared 
regulator, as was mentioned. 

Consistency with neglecting the LPM contribution re¬ 
quires to take the large JVc limit. The production rate is 
then given by Eq. (4.4) with the replacement — 1 —>■ 
in the prefactor. By extrapolating this expression to 
the Nc = 3 case, and inserting In from Eq. (2.11), we 
have 


E 


dr 

d^p 


Comp 


47r4 


(48) 

/ 


with the dimensionless function 


-[1-1-2 cos(/3n(5)] 


2ET 


n—l 
cco 


pCO 

/ dsG{s,m})\ 

Jo 


(4.9) 


— nl3{m^E/s+S/ (4i5)) 


B. Pair annihilation 


In a way similar to the case of the Compton scattering 
process, the photon production rate due to pair annihi¬ 
lation can be written as: 


E 


dr 

d^p 


2{2tt) 


— T 


f d^k 

CO 

f d^p' 

J (27r)3j 

' (2^)3 j 

' (27r)3 


|M1 = 


yj-(27r)"^i5*'"’‘^(fc + k' — p — p') 

Ek') [1 + nBcdi 
r°° r 


SEkEk'Jp' 

X nFa{E\s) nFh(£’k') [1 + nBcd{\p' 
-1 T 


o-PE 


32(27r)6 E 
X g*/3(Qa-Q6-Qc-|-Qd) 

X In 1 — g“/5'*/(4^^)g*^(Q=-Qd) 


‘+ ~\M} 

dt -—^ 


(4.10) 
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where we have used again the Boltzmann approximation 
in the initial state and introduced u = (jp — k'Y, u = 
u — and i± = —s[l =F j^j^. The squared 

matrix element summed over spins is: 




p {tu — 2nYi — AmY 


+ i ~ 2m^M — (s — 4m^) 

tu 


(4.11) 


Using s + t + £t = 0, the result reads 


E 


dr 

#p 


aUs T 


-PE 


J2q} ds 


4m2 


oo ^ 

I ll - 1) H{s, nY), 

n—1 

(4.12) 



FIG. 8: The quantities X and Y as functions of T with E — 2 
GeV. The parameter set A is used. Shown for comparison are 
results obtained in the leading-log approximation (LL) and in 
the completely deconfined phase (Adeconf, Ldeconf)- 


where 


H{s,m?) = ( 1 -h 


4m^ 


+ 1 


\ 1 4m2 

Inserting Nc = 


c.dr 

00 

(N 



X In 


1 -h — Aw? js 


1 — y^l — Am? I s 


(4.13) 


In the deconfined phase with In = 1, Eq. (4.12) reduces 
to 


pair j 


with the dimensionless function 


y(T) = 


AET 


4m2 


dsj2- 

n—1 


^-snP/(4E) 


X [2cos{n/3Q) -h 1]^ E[{s,m^). 


(4.17) 


(4.18) 


E 


dr 

#p 


aa<, T 


pair 


167r4 E 




E9/(^c-1) 


[ ds H{s,ttY). 

J 4m2 ^ ' 

(4.14) 


In the limit —>■ 0, we have 

dr 


E 


d^p 


aug 

Att^ 

X In 




n—1 


ET 

2" ■ 

m ^ 


(4.15) 


which again agrees with the leading-log result [17]. 
In the large limit, Eq. (4.12) becomes 


E 


dr 

#p 


attg T 


-PE 


167r4 E 




ds 


4m2 


X V - e-'<'^/5/(4iS)^2 ;2 ^2) ^ 


(4.16) 


C. Numerical results 

We proceed by first analyzing separately the con¬ 
tributions from Compton scattering and pair anni¬ 
hilation. In Eqs. (4.8) and (4.17), the prefactors 
aasT'^e~^^Y^jq^f/{A-KY are common, so we plot the T- 
dependent functions X and Y in Fig. 8. The energy is 
set to E = 2 GeV, and the parameter set A is used to 
generate m(T). As Polyakov loop input a fitting func¬ 
tion for the lattice data of 1{T), plotted in Fig. 1, has 
been used. Comparing these quantities with the corre¬ 
sponding ones at the leading-log order (Eqs. (4.7) and 
(4.15) at large Nc limit), we see that the T dependence 
can be well described by the leading-log approximation 
when T is larger than about 250 MeV. At such tempera¬ 
tures m becomes negligible which makes the leading-log 
approximation more accurate. Also shown in the same 
figure are X and Y in the deconfined limit with I = 1 
(Eqs. (4.6) and (4.14) at large Nc limit). One observes 
that the Polyakov loop effect suppresses these quantities 
significantly at all T. Note that the contribution from 
pair annihilation is larger than that from the Compton 
scattering process in the limit 1 = 1, whereas this is not 
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FIG. 9: The functions X + Y and Z at various temperatures T 
with E = 2 GeV. The parameter set A is used. The solid (red) 
line represents X + Y while the dotted (blue) line represents 

Z. 


FIG. 10: The sum X + Y calculated using the parameter sets 
A, B, C,D as a function of T at i? = 2 GeV. 


the case when incorporating the Polyakov loop. This be¬ 
haviour derives from the fact that X is linear in In while 
Y involves (/n)^, which makes the suppression of Y due 
to the Polyakov loop more significant. 

In order to examine the effect of the dynamical quark 
mass we compare results with finite m to those obtained 
using m = 0. With inclusion of m{T) the total photon 
production rate is given by 

= + ( 4 . 19 ) 



In the m = 0 case the infrared singularity is regularized 
by the thermal quark mass. The calculation requires a 
resummation, the derivation of which is given in the ap¬ 
pendix. Here we just quote the result, adding Eqs. (A.7), 
(A.8), and (A.23): 




massless 




/ 


(4.20) 


where 


Z{T) = TT^ - 89 -f 


9 


In 




9 


- TT^ 159^ - 89 + r: - 2.5087 


00 1 

— 2 E (3(—-I- 2) cos(27rn9) -|- cos(47rn9) 


_ .69 + 692 )ln(---9+-9^ 


+ ( - - 29 ) In 


1 _ 9 

16 2 


(4.21) 


which is found using the expression for large N,. and in¬ 
serting Nc = 3. Here 9 = Q/{2ttT). Note that although 
an arbitrary cutoff /i is introduced in the derivation, the 
result is independent of /i. 

The T-dependent functions X + Y and Z are shown in 
Fig. 9. One observes that A + F is significantly smaller 
(larger) than Z when the temperature is smaller (larger) 
than about 190 MeV, close to the transition temperature 
To of parameter set A. This behavior of A + F can be 
understood as follows. In the low temperature region, m 
is of order or larger than T. From the exponential part of 
Fq. (4.4) it follows that the dominant contribution to A 
comes from the region s ~ Em so that the suppression 
factor results in A. As for the pair annihilation 

contribution, the lower bound of the integral in Fq. (4.12) 
gives the exponential suppression factor e~^ /{et) y 
when m is large. On the other hand, when T exceeds 
about 190 MeV, m becomes small and makes the leading- 
log approximation (Fqs. (4.7) and (4.15)) more accurate. 
This leads to an enhancement of A + F compared to Z 
because of the factor \w{ET/m?). 

While these considerations remain valid, the detailed 
high-T behavior needs nonetheless to be examined with 
greater care. When m is much smaller than the thermal 
quark mass, which is realized for T > 200 MeV according 
to Fig. 2, this thermal mass works as an infrared regula¬ 
tor instead of m, and thus the role of the leading factor 
\iv{ET/m?) appearing in A + F is now to be replaced by 
ln(T/(9^T)) in Z. At temperatures well above 200 MeV 
it is therefore appropriate to use Z rather than A + F. 
In summary, for practical applications we suggest to use 
A + F for T < Tm and Z for T > T^, where T^ is the 
temperature at which A + F equals Z. 

Next we examine how variations of the temperature Tq 
and of the steepness 6 of the chiral crossover transition 
influence the photon production rate. Fig. 10 shows the 
quantity A+F, using the parameter sets A, B, C, and D, 
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as functions of T. We focus on the low T region (T < T^) 
where the results including the dynamical quark mass 
m are expected to be reliable. Varying S changes the 
photon production rate only moderately, as can be seen 
by looking at X + Y using the parameter sets A, C, and 
D with Tq = 180 MeV. By contrast, varying Tq induces 
a more significant change in V + F as can be seen by 
using the parameter set B with Tq = 150 MeV. The 
shift of the quark mass profile m{T) to a lower transition 
temperature implies an increase ot X + Y already around 
T Tq where the quark mass starts to drop rapidly. 

V. SUMMARY AND CONCLUDING REMARKS 

In this work dilepton and photon production rates in 
the quark-gluon phase close to the transition tempera¬ 
ture Tc have been calculated using a model that takes 
into account effects of confinement and spontaneous chi¬ 
ral symmetry breaking (ySB), mainly focusing on the 
latter. Implications of ySB are realised in the form of 
dynamically generated quark masses as the temperature 
approaches Tc from above. Through its induced kine¬ 
matic constraint, this effect significantly suppresses the 
production of dileptons with large three-momenta. For 
dileptons emitted back-to-back with zero momentum the 
production rate is reduced more modestly. 

We have also examined the dependence of the dilepton 
production rate on the parameters that characterise the 
chiral crossover transition, suggesting that information 
about the transition can indeed be encoded in the pro¬ 
duction rate of dileptons with large three-momentum. In 
particular variations of the transition temperature imply 
different spatial momenta at which the dilepton produc¬ 
tion vanishes as a consequence of phase space suppres¬ 
sion. 

A comparison of our calculated dilepton production 
rate with the results of other nonperturbative calcula¬ 
tions is useful at this point. The dilepton rate (Eq. (3.5)) 
for V/ = 2 at IpI = 0 and T = 240 MeV is plotted as a 
function of in Fig. 11. The parameter set A was used 
to calculate m(T). Note that in this casep° = M = \/^, 
the invariant mass of the dilepton, and we use M as the 
variable in this hgure. For comparison the results from a 
PNJL model calculation [33], lattice QCD [34], and a per¬ 
turbative calculation are also shown. The perturbative 
(Born approximation) result is obtained using Eq. (3.5) 
with A = B = C = 1. We note that our expression 
for the dilepton production rate (Eq. (3.7)) has the same 
functional form in terms of I and m as that^ obtained by 
using the PNJL model [33]. This is not surprising since 


From Eqs. (4.36) and (4.46) in Ref. [33], one can reproduce 
Eq. (3.7) in our paper by considering the case of zero chemical 
potential. To confirm this, one needs to rewrite f{Ep) — 1 in 
the two equations in Ref. [33] to —f{Ep — lu) by changing the 
integration variable p. 


the leading dilepton production amplitude {qq —>■ 7 —>■ (I) 
is the same in both approaches. Beyond leading order, 
loop corrections enter, which make the results depend on 
details of the model used. 

At T = 240 MeV the dynamical quark mass rn{T), in 
our approach as well as in the PNJL model, is negligi¬ 
ble compared with T. The difference in Fig. 11 between 
our result and the PNJL model comes solely from the 
difference of the Polyakov loop values used in these cal¬ 
culations. Our smaller Polyakov loop stems from sub¬ 
tracting the perturbative correction, so that our dilep¬ 
ton production rate is enhanced compared to that of the 
PNJL model, as explained in Ref. [9, 17]. The overall 
enhancement of the dilepton production rate induced by 
the Polyakov loop is also explicit in the fact that the per¬ 
turbative (Born) rate is smaller than both the PNJL and 
semi-QGP results. 

At small invariant mass M one observes that our cal¬ 
culated dilepton rate is significantly lower than the one 
deduced from lattice QCD. One should note that the lat¬ 
tice results display considerable sensitivity with respect 
to the ansatz chosen for the form of the vector current 
spectral function in the low-mass region, as discussed in 
Ref. [33]: a “free-field” ansatz [35] implies a very small 
dilepton production rate at small M, whereas a “Breit- 
Wigner” ansatz [34] yields a large dilepton rate reflecting 
substantial contributions from hadronic sources. 

The photon production rate is found to be sensitive 
to ySB effects when T is less than about 190 MeV. We 
have studied the parameter dependence of this rate. It is 
found to vary significantly at low T around T^ depending 
on the choice of the transition temperature. 

Comparing our calculated photon rate with corre¬ 
sponding results using other methods is once again of 
some interest. Fig. 12 shows our result, Eq. (4.19), as 
a function of the photon energy E in comparison with 
the semi-QCP photon rate [9, 17] taken at leading-log 
order and using m = 0, and with the perturbative QCP 
result [28, 29]. Here we have set T = 180 MeV, V/ = 2, 
and used the parameter set A. The large-Ac limit is taken 
consistently in all three case studies. The semi-QCP ap¬ 
proach at leading-log order and m = 0 corresponds to 
Eq. (4.20) with 

Z = (5.1) 

while the perturbative QCP result is given with Z in 
Eq. (4.21) at g = 0. As seen in Fig. 12 the perturbative 
photon rate is larger than the semi-QCP rate over the 
whole energy range. The reason is that both Polyakov 
loop and constituent quark mass effects suppress the pho¬ 
ton production as discussed in Sec. IV C and Ref. [9, 17]. 
At the lower energies (T ~ 1 — 2 CeV) our result ex¬ 
ceeds the one calculated at leading-log order and with 
m = 0. This difference comes primarily from the term 
beyond leading-log. At high photon energies the leading- 
log form is approached asymptotically. Indeed, at suf¬ 
ficiently large E, the expressions for X and Y can be 
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approximated by their leading-log forms, Eqs. (4.7) and 
(4.15), so that X + Y equals Z as given in Eq. (5.1), ex¬ 
cept for the replacement ln{E/{g^T)) —>■ hi{ET/m'^). At 
large E, this difference is negligible, and so the photon 
production rates in the two approaches become approxi¬ 
mately equal. 

A few concluding remarks are in order concerning 
possible ySB effects in high-energy heavy-ion collisions 
(HIC). Dilepton production in HIC is commonly dis¬ 
cussed using the invariant mass M instead of |p|. With 
inclusion of dynamical quark masses, the kinematical 
constraint for dilepton production reads M > 2m{T). 
In the NJL model the dynamical quark mass is still quite 
large (m ~ 300 MeV) at temperatures around T ~ 150 
MeV, so dileptons with M < 600 MeV are forbidden in 
this temperature range. At slightly higher temperatures 
where the dynamical quark mass decreases the kinematic 
constraint is correspondingly less stringent. While such 
effects can in principle be observed in HIC, one recalls, 
on the other hand, that the low-mass dilepton spectrum 
with M < 0.8 GeV results primarily from hadronic cur¬ 
rents involving pseudoscalar and vector mesons rather 
than liberated quark-antiquark pairs [36]. In this range 
of invariant masses, dilepton spectra are governed by 
hadronic sources. 

An interesting point of discussion is the photon ellip¬ 
tic flow V 2 observed at RHIC [10]. As demonstrated in 
Ref. [9], the suppression of the photon production rate 
from the QGP phase due to the confinement effect can 
increase the total elliptic flow, a welcome feature in order 
to understand the unexpectedly large V 2 for photons [9]. 
Our results suggest that the effect of ySB can add to 
increasing the total photon V 2 further at T between 150 
and 190 MeV by its mechanism of suppressing the pho¬ 
ton production in that temperature range. Of course, 
other sources of photon production in HIG also need to 
be considered, such as photons coming from the initial 
state [14], from the thermalization process [15], and from 
the hadron phase [11]. A complete evaluation of the 
photon V 2 requires summing the contributions from all 
the photon sources in extended hydrodynamic simula¬ 
tions [9, 12, 13]. 
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Appendix: Photon Prodnction Rate when m = 0 

In this appendix we derive the photon production be¬ 
yond leading-log order in the limit of vanishing quark 
mass, TO = 0. The result at leading-log order has been 
reported previously [9, 17], but the derivation beyond 
this order has not been performed yet. The infrared sin- 



M [GeV] 

FIG. 11: The dilepton production rate for Nf = 2 at jpj = 0 
and T = 240 MeV as a function of the dilepton invariant 
mass M. The solid (red) line represents our result (Eq. (3.5)), 
the dashed (green) line is the result of a PNJL model cal¬ 
culation [33], the dotted (blue) line is the result from lat¬ 
tice QCD [34]. The dash-dotted (magenta) curve shows 
the perturbative (Born) dilepton rate (using Eq. (3.5) with 
A^B^C =1) . 



E [GeV] 

FIG. 12: The photon production rate for Nf — 2 sk, T = 180 
MeV as a function of the photon energy E. Solid (red) 
line: present result (Eq. (4.19)); dashed (black) line: semi- 
QGP model with m = 0 and using leading-log approxima¬ 
tion [9, 17]; dotted (blue) line (labeled QGP with m = 0 and 
Q = 0): perturbative QCD calculation including resumma¬ 
tions [28, 29]. 


gularities in the contributions from Gompton scattering 
and pair annihilation require introducing a cutoff /r as 
a regulator. In principle g is an arbitrary quantity, but 
it should satisfy gT <C /r T. The phase space inte¬ 
gral appearing in the Gompton scattering process and 
the pair annihilation contributions are restricted by the 
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FIG. 13: The photon self-energy at the one-loop level. The 
solid line with a blob represents the HTL-resummed quark 
propagator. 


cutoff /r [29]: 2^^ < s and + < t < In the 

phase space excluded by /i, the momentum of the ex¬ 
changed quark is smaller than /r, and thus much smaller 
than T. Such soft quarks should be treated using the 
hard thermal loop (HTL) resummed propagator [19, 37], 
which contains the information of the thermal mass of 
the quark, so we calculate the resummed one-loop di¬ 
agram drawn in Fig. 13 in order to obtain the photon 
production rate due to the soft momentum exchange. 

The hard parts from Compton scattering and pair an¬ 
nihilation are given by Eqs. (4.1) and (4.10), setting 
TO = 0, and the phase space integration is constrained 
by the cutoff /i as explained above. The resulting expres¬ 
sions are as follows [17]: 
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Here we have used the Boltzmann approximation for the 
initial state. By performing these integrations, we get 
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In the large-A^c limit these expressions become 
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With Nc = 3 they reduce to: 
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where Cciaisher — 1.2824 is the Glaisher constant. 
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On the other hand, the contribution from the soft part 
reads [17] 



X [nFa(w - E) - npaioj)] 

X Tr [p*(w, k)p(a; - k - p)], 

(A.13) 


which corresponds to the one-loop diagram drawn in 
Fig. 13. The summation includes the color indices a. 
Here the bare (HTL-resummed) quark spectral functions 
are: 


p{u}, k) = 27re(w) lf:5{k^) 
7° — 7 • k 


Pl(w,k) = 




7 • k 


(A.14) 

Pa-(w,k), 

(A.15) 


where e{uj) is the sign function. The HTL spectral func¬ 
tion after decomposition reads 

p*±(w,k) = 27r[Z±a(|k|)5(w - a;±a(|k|)) 

+ ^=Fa(|k|)<^(w -b W;pa(|k|))] (A.16) 

Here Z±a[\^\) is the residue, a;±Q(|k|) the dispersion re¬ 
lation, and ma(T) is the thermal quark mass. Their re¬ 
spective expressions are: 


where Qo(a:) = (ln(a;-|-l)/(a;—1))/2. The explicit form of 

spacelike • • i - • i • j ^ -j. 

PcL± irrelevant in our analysis, so we do not write 

it here. The thermal quark mass is given by [19] 




4 

1 



Nc 

qb\ + Nc 

. 6=1 




(A.19) 


where qa = Qali^irT). At large Nc, the last term is ne¬ 
glected. The production rate is reduced to the following 
form [17]: 


E 


dr 

d^p 


4a 


soft 


(27r)3 


,-/3B 




-l + ln^ + Jy/ d\k\ 


mi: m 


a 


X W+adk]) - W-a(|kl) - 


(A.20) 


The integration in the equation above does not depend on 
a, and its numerical value approximately —0.16258 [28]. 

Using Eq. (A.19) in the large-A^c limit, and considering 
the case of A^c = 3, we get 


Z±a{m) = 

w±a(|kl)T|kl 


^L(|kl)-|kl = 


2ma 


(A.17) 


ml 


1 ^ oUeTbMl ± 1 


(A.18) 


9 9 , 9 . 

ml = ml = (1 - 4g -b 4q^) , 


ml = 


16 

Zg'^T'^ 

16 


1-3,1. 


With these expressions one finds: 


(A.21) 

(A.22) 
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dr 


soft 


qf— 

47r4 2 

/ 



1.32516) I 




(A.23) 


The final result in the large limit is given by the sum Eqs. (A.7), (A.8), and (A.23), when = 3. 
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